Localization in Correlated Bi-Layer Structures: 
From Photonic Cristals to Metamaterials and Electron Superlattices 
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In a unified approach, we study the transport properties of periodic-on-average bi-layered photonic 
crystals, metamaterials and electron superlattices. Our consideration is based on the analytical 
expression for the localization length derived for the case of weakly fluctuating widths of layers, 
that also takes into account possible correlations in disorder. We analyze how the correlations lead 
to anomalous properties of transport. In particular, we show that for quarter stack layered media 
specific correlations can result in a tj 2 -dependence of the Lyapunov exponent in all spectral bands. 

PACS numbers: 42.25.Dd, 42.70.Qs, 72.15.Rn 
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Introduction. In recent years much attention was 
paid to the propagation of waves (electrons) in periodic 
one-dimensional structures (see, e.g. [l| and references 
therein). The interest to this subject is due to vari- 
ous applications in which one needs to create materials, 
metamaterials or electron superlattices with given trans- 
mission properties. One of the important problems that 
still remains open, is the influence of a disorder that can 
not be avoided in experimental devices. Such a disor- 
der can be manifested by fluctuations of the width of 
layers and distance between layers, or due to variations 
of the medium parameters, such as dielectric constant, 
magnetic permeability or barrier hight (for electrons) 

As is well known, the main quantity that absorbs the 
influence of a disorder is the localization length ent irely 
determining transport properties in a ID geometry [lOj. 
In contrast to many studies of the wave (electron) propa- 
gation through random structures, mainly based on var- 
ious numerical methods, in this Letter we develop an 
analytical approach allowing us to derive the unique ex- 
pression for the localization length that is valid for pho- 
tonic crystals, metamaterials and electron superlattices. 
Another key point of consideration is that we explicitly 
take into account possible correlations within a disorder, 
that may be imposed experimentall y. As was recently 
shown, both theoretically 11, 12, 13, and experimen- 
tally pi, [lj|, specific long-range correlations can signifi- 
cantly enhance or suppress the localization length in de- 
sired windows of frequency of incident waves. 

Model. We consider the propagation of electromag- 
netic wave of the frequency u> through an infinite array 
of two alternating a and b— layers (slabs). The slabs are 
specified by the dielectric constant e a ,b, magnetic perme- 
ability /V&, refractive index n a ,& = Je a ,bp>a,bi impedance 
Z a ,b = \Z^a,b/s a ,b and wave number k a .b = um a ,b/c. We 
assume that the z-axis is directed along the array of bi- 
layers perpendicular to the stratification. Within the lay- 



ers, the electric field obeys the wave equation, 
d 2 



dz 2 



+ k a,b ^a.b{z) = 0, 



(1) 



with two boundary conditions on the interfaces z = Z{ be- 
tween slabs, ipa(zi) = ipb{zi) and n~ 1 vb' a (z i ) = n^ip'^zi). 

A disorder is incorporated in the structure via the 
random widths of the slabs, a(n) = a + g a (n), b(n) = 
b + Qb(n). Here n enumerates the elementary a6-cells, a 
and b are the average widths of layers and g a {n), Qb{n) 
stand for small variations of the widths. In the absence 
of disorder the array of slabs is periodic with the pe- 
riod d = a + b. The random sequences g a ,b{n) are sup- 
posed to be statistically homogeneous with zero average, 
{Sa,b{n)) = 0, and the correlations are fully determined 
by the binary correlation functions 

{ ej {n) ej {n')) = {g 2 {n))K 3 {n - n'), j = a,b (2) 

(Qa(n)Qb(n')) = (g a (n)g b (n))K ab {n - n ). (3) 

In what follows the average (...) is performed over the 
whole array of layers or due to the ensemble averaging, 
that is assumed to be the same. The two-point auto- 
correlators K a ^{n — n') as well as the inter-correlator 
K a b(n — n') are normalized to one, K a ^(0) = K a b(0) = 
1. The variances {g 2 ab (n)) are of positive value, while 
(g a (n)gb(n)) can be both positive and negative. Note 
that \(g a {n)gb(n))\ = y/ (g 2 (n)) {g 2 {n)). We assume the 
positional disorder be weak, k 2 b (g 2 lb (n)) <C 1, allowing 
us to use an appropriate perturbation theory. In this 
case, all transport properties are entirely determined by 
the randomness power spectra K, a {k), K.b(k), and K, a b(k), 
defined by the relation, fC(k) — J27L-oo K{r) exp(— ikr). 
All the correlators K(r) are real and even functions of the 
difference r = n—n'. Because of this fact and due to their 
positive normalization, the corresponding Fourier trans- 
forms IC(k) are real, even and non- negative functions of 
the dimensionless lengthwise wave- number k. 
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Method. Our aim is to derive the localization length 
(LL) loo(u>) in the general case of either white or colored 
disorder. On the scale of individual slabs the solution of 
Eq. ([1]) can be presented in the form of two maps for nth a 
and b layer, respectively, with corresponding phase shifts 
<?aA n ) = ¥>a,b + €a,b( n )> where (p a = k a a, ip b = k b b, and 
^a,b(n) = k a ,bQa,b(n)- By combining these maps with the 
use of the boundary conditions, one can write the map 
for the whole nth elementary afr-cell, 



Vn+i = -C n x n + D n y n . (4) 



Here x n = tp a (zan) and y n = k~ l ij)' a {z an ), the index n 
corresponds to the left edge while n + l^stands for the 
right edge of the nth cell. The constants A n , B n , C n , D n 
depend on tp a (n), <pb{n) and on Z a /Z b = k b (i a /k a (i b . 

Eq. ([4]) can be treated as the map of a linear oscillator 
with time-dependent parametric force [ijj • Without dis- 
order the trajectory x n , y n creates an ellipse in the phase 
space (x,y), that is an image of the unperturbed mo- 
tion. It is convenient to make the transformation, x n = 
v~ 1 Qn cos r — vP n sinr, y n — v~ 1 Q n sin T + vP n cos r, to 
new coordinates Q n ,P n , in which the unperturbed tra- 
jectory occupies the circle, Q n +i = Q n cosj + P n smj, 
P n +i = —Qn sin 7 + P n cos 7, in the phase space (Q,P). 
Here r and v can be found from Eq. ([4]), and 7 determines 
the Bloch wave number k = "f/d arising in the relation 
ip(z + d) = cxp(iKd)ip(z) for the periodic array, 

cos 7 = cos (p a cos (p b ^- + ^r) sin sin 93;,. (5) 
I \Z a Z b J 

In order to take into account the disorder, we expand 
the constants A n , B n , C n , D n up to the second order in 
the perturbation parameters £ a (n) <C 1 and £&(n) <C 1. 
Then, one can transform the variables x n , y n into Q n , P n . 
After getting the perturbed map for Q n ,P n , we pass to 
action-angle variables R n , 9 n via the standard transfor- 
mations, Q n = Rn cos 6 n , P n — RnSinOn- This allows 
us to derive the relation between R n +i and R n keeping 
linear and quadratic terms in the perturbation, 



R 



■71+ 1 



1 + Un)Va(n) + Zb(n)V b (n) - C Q 2 (n) - 4 2 (n) 



+e a (n)W a (n)+e b (n)W b (n)+Un)Un)W ab (n), (6) 

where V a (n), V b (n), W a (n), W b (n), W ab (n) are compli- 
cated functions of 9 n and the model parameters. 

Localization length. The LL can be expressed via the 
Lyapunov exponent (LE) A = d/loo{u>) defined by 17] . 



\ 1 /1 ( 



(7) 



In deriving the LE it was assumed that the distribution of 
n is homogenous within the first order of approximation. 
This assumption is correct apart from the band edges, 



7 = 0;7r, and the vicinity of the center, 7 = it/2 [13j |. 
The calculation of the LE has been done with the use of 
the method developed in Refs. [TllElEi. 

Omitting all 

details, here we refer to the final result for the total LE, 



A = - m „ cr 2 /C a (27)sin 2 <y9 6 + cr fc 2 /C b (27)sin 2 ( ( 9 a 
8 sin 7 L 

-2a 2 b /Ca6(27)sin^ a sin^ 6 cos7 , (8) 



k a k b (g a (n)g b (n)), 



where <6 = k ib(e 2 a,b( n ^\ a lb 
and vj = Z a /Z b — Z b /Z a is the mismatching factor. 

The expression ((8]) generalizes the results obtained in 
Refs. [H, 0, HH for particular cases, and is in a com- 
plete correspondence with them. Let us now discuss the 
derived expression in some applications. 

Conventional Photonic Layered Media. In this case all 
the optical characteristics, e a ,b, l^a,bi n a ^ b , Z a j, are pos- 
itive constants. One can see that if the impedances of 
a, 6-slabs are equal, Z a = Z b , the mismatching factor w 
entering Eq. ||5J) vanishes and the perfect transparency 
emerges (A = 0) even in the presence of a disorder. Since 
the layers are perfectly matched, this conclusion is gen- 
eral [l|, and does not depend on the strength of disorder. 
According to Eq. (JSJ) , in such a case the stack-structure is 
effectively equivalent to the homogeneous medium with 
the linear spectrum, n = cun/c, that has no gaps, and 
where the refractive index is n = (n a a + n b b)/ (a + b). 

For the Fabry-Perot resonances appearing when uj/c = 
s a ^/n a a and oj/c = s b ir/n b b, with s a , b = 1, 2, 3, ... , the 
factors siniy9 a and sin^b in Eq. (jSJ vanish, thus giving 
rise to the resonance increase of the LL. In a special case 
when n a a/n b b — s a /s b some resonances from different 
layers coincide and would give rise to the divergence of 
the LL. However, such a situation can arise only at the 
edges of spectral bands, where 7 = 0, n and sin 2 7 also 
vanishes. Therefore, the LL gets a finite value at these 
points instead of diverging (due to resonances) or van- 
ishing (due to band edges). The above statement is not 
valid at the bottom of the spectrum, u>/c= 0, where the 
analysis has to be done separately. Note, nevertheless, 
that for a white noise, /C a ,b = 1, fcab = 0, ±1, the LE 
obeys the conventional dependence Atxw 2 for u> — > 0. 

A special interest is in the long-range correlations lead- 
ing to the divergence of the LL in the controlled win- 
dows of frequency u>. This effect is similar to that found 
in more simple ID models with the correlated disorder 
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lore simple 



In our model this effect is due to 



possibility to have the vanishing values of all Fourier 
transforms, K a ,b = fcab = 0, in some range of frequency 
lu. This fact is important in view of experimental re- 
alizations of random disorder with specific correlations. 
In particular, one can artificially construct an array of 
random bi-layers with such power spectra that abruptly 
vanishes within prescribed intervals of w, resulting in the 
divergence of the LL. Also, specific correlations [19( in a 
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disorder can be used to "kill" a sharp frequency depen- 
dence associated with the term sin 2 7 in the denominator 
of Eq. ([5]). It is noteworthy that in the middle of spectral 
Bloch bands (7 = 7r/2), the third term vanishes and the 
inter-correlations do not contribute to the LL. 

The typical dependence A(w) for the conventional pho- 
tonic bi-layer stack is shown in Fig. 1. 



within any pair of a, 6-layers, although the width of both 
a, 6-layers fluctuates randomly. For such correlations, the 
quadratic w-dependence seems to remain within the non- 
perturbative regime as well. 





FIG. 2: Lyapunov exponent versus frequency for quarter 
stack layered medium, for the "plus correlations" (left) and 
"minus correlations" (right), see Eqs. |TT} and (|9}. Here 
FIG. 1: Lyapunov exponent A versus frequency w in arbitrary (Z a + Z b ) 2 /Z a Z b = 2.4, Z(Q 2 a (n))(irn a / c) 2 = 3.2, Ka = 1. 
units for vj 2 n 2 a {gl{n))/{2c 2 ) » 12.28, vo 2 n 2 b ($j(n))/(2c 2 ) « 
0.27 and n a a/c — 1.6, n b b/c = 0.4. Left: photonic layered 
medium, Eqs. |H} and JSJ. Right: the RH-LH bi-layers with 
"plus" in Eq. 



Quarter Stack Layered Medium. This term is typically 
used when two basic layers, a and b, have the same optical 
width, n a a — n^b (see, e.g. Since ip a = ifb, in this 

case the dispersion relation ((SJ) takes the form, 



cos 7 = 1 

' 2 



1 (Z a + Z h f 



Z a Zb 



sin 2 (k a a). 



(9) 



One can see that starting from the second band the top 
of every even band coincides with the bottom of next odd 
band at 7 = 0. The gaps arise only at 7 = n. 

With the use of Eq. (jQj) the LE can be written as 



Metamaterials. A special interest is in the mixed sys- 
tem in which the a-layer is a conventional right-handed 
(RH) material and &-layer is a left-handed (LH) material. 
This means that e a ,/i a ,n a > 0, whereas Eb,[tb,nb < 0. 
However, both impedances remain positive, Z a ,Zb > 0. 
Remarkably, in comparison with the conventional stack- 
structure, the expression (JSJ) for the LE in this case stays 
the same. The only difference is that in the dispersion 
equation J5} the sign "plus" has to substitute for "minus" 
at the second term (the phase ipb = kbb = —uj\rib\b/c). 
Such a "minor" correction can drastically change the fre- 
quency dependence of LE, see Fig. 1. Nevertheless, the 
LE caused by positional disorder, typically obeys the con- 



ventional dependence, A 



when 1 



0. 



A = 



2/(7) 



Z = 



{Z a ~ Zbf 



(10) 



8 cos 2 ( 7 /2)' ~ Z a Z b 
/( 7 ) = CT X(27) + <#C 6 (2 7 )- 2^ 6 /C a6 (2 7 ) cos 7. 

Thus, the LE is finite or vanishes at 7 = and diverges 
at 7 = it. 

It is instructive to analyze the simplest case of corre- 
lations when either £ a (n) = £b(n) (plus- correlations) or 
£a(n) = —£,b(n) (minus- correlations). In this case one 
gets /( 7 ) = 2cr 2 /C a (27)(l =F COS7), correspondingly. As a 
result, for the LE one can obtain, 



A_ 



Za 2 JC a (2j) 



(11) 



As one can see, A + oc lu 4 at the bottom of the spectrum 
(u) — ► 0), in contrast to the conventional dependence A oc 
lu 2 . Another non-conventional dependence, A cx to was 
recently found Q in a different layered model with left- 
handed material. It is interesting that for the minus- 
correlations and /C a (2 7 ) = 1, the w-dependence of the LE 
is quadratic for any energy inside the spectral bands (see 
Fig. 2). In this case the total optical length is constant 



Note also that the ideal mixed stack (e a = fx a = n a = 
1, Eb = Hb = Tib = — 1, Z a = Zb = 1) has perfect trans- 
mission, A = 0, independently of a positional disorder. 

One of the interesting features of the mixed layered 
structures is that for the RH-LH quarter stack (n a a = 
\nb\b) the average refractive index vanishes, n = (n a a — 
\nb\b)/(a + b) — 0. As a consequence, it follows from 
Eq. ([5]) that the spectral bands disappear and therefore, 
the transmission is absent, apart from a discrete set of 
frequencies where k a a — sir and 7 = (s = 0, 1, 2, . . .). 
Evidently Eq. ([5]) is not valid in such a situation and an 
additional analysis has to be done. 

It is important that in reality there is a frequency 
dependence of the permittivity Eb(ui) and permeability 
/ifc(w) [1]. This fact is crucial in applications. In par- 
ticular, it leads to the following peculiarities. First, the 
mismatching factor w in Eq. ([8]) can vanish for specific 
values of frequency u> only, thus resulting in a resonance- 
like dependence for the transmission. Second, for typi- 
cal frequency dependencies the refractive index of 6-slabs 
takes an imaginary value giving rise to the emergence 
of new gaps. Specifically, such a gap can arise at the 
origin of spectrum, lu = 0, in contrast to conventional 
photonic crystals. It can be seen that in many aspects 
the wave transport through the bi-layered metamaterials 



4 



is resembling to that of the electrons through double- 
barrier structures. 

Electrons. The developed approach can be also ap- 
plied to the propagation of electrons through the bi- 
layer structures with alternating potential barriers of the 
amplitudes U a and Ub and slightly perturbed widths. 
Indeed, the stationary ID Schrodingcr equation for an 
electron with effective masses m a , wi& inside the barri- 
ers and total energy E can be written in the form of 
Eq. ([l]), in which the partial wave numbers are asso- 
ciated with the barriers, k a = yj2m a {E — U a )/h and 
kb = w2mb(E — Ub)/h. Another change, n a b — * m a ,6, 
should be done in the boundary condition on hetero- 
interfaces, m^ 1 ij)' a {zi) — m'^ 1 ijj' b (zi). Correspondingly, 
Zb/Z a = k a mb/kbm a in the dispersion relation |(5J) and 
in the expression © for the LE. 

If the energy E is smaller then the heights of both bar- 
riers, E < U a , Ub, the electron wave numbers are purely 
imaginary. As a consequence, the electron states are 
strongly localized and the structure is non-transparent. 

For other case when U a < E < Ub, the tunneling propa- 
gation of electrons emerges. In this case the wave number 
k a is real while kb is imaginary. Therefore, the electron 
moves freely within any a-barrier and tunnels through the 
^-barriers. Thus, the expressions ijHJ) for A and (O for 7 
have to be modified according to the change, kb — > i\kb\ 
and sin(kbb) — ► i sinh(|fci,|6). As a result, the Fabry-Perot 
resonance increase of the LL arises only due to the second 
and third terms of Eq. ([8]). The general expression for 
the LE can be essentially simplified for a particular case 
of an array with delta-like potential barriers and slightly 
disordered distance between them. The corresponding 
expression for the LE is in full correspondence with that 
obtained in Ref. 1121 1. 
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FIG. 3: Lyapunov exponent vs. energy for electrons in bi- 
layered structure. Here 2m a /~h 2 = 2mt//i 2 = 1, U a — 0.15, 
U b = 1.6, a = 0.35tt, b = 0.65tt, (gl(n)} = (0.1a) 2 , 
(f?!>( n )} — (0.2&) 2 The transition between tunneling and over- 
barrier scattering occurs at E/n = 0.5. 

For the over-barrier scattering, when U a < Ub < E, 
both wave numbers, k a and kb, are positive and the elec- 
tron transport is similar to that for the conventional pho- 
tonic stack but with dispersive parameters. The example 
of the energy dependence of the LE is given in Fig. 3. It 
is interesting that if E — (Ubm a — U a mb) / (m a — mj), 
an electron does not change its velocity in the barri- 
ers, hk a /m a = hkb/nib, although its momentum changes, 



hk a > hkb- The LE vanishes in this case. Thus, such a 
"free" electron motion is equivalent to that in a homoge- 
neous medium with perfect transmission. 

Conclusion. We derived the expression for the inverse 
localization length for quasi-periodic bi-layer structures, 
whose widths are weakly perturbed. Our result can be 
applied both to conventional photonic crystals and to 
metamaterials, as well as to the electron superlattices. 
Another feature of the approach is that it takes into ac- 
count possible correlations in a disorder that can lead to 
anomalous frequency (energy) dependence of transport 
properties. Due to the correlations one can significantly 
enhance or suppress the transmission/reflection through 
the bi-layered devices within the prescribed windows of 
frequency (energy) of electromagnetic (electron) waves. 
Our results may have a strong impact for the fabrica- 
tion of a new class of disordered optic crystals, left /right 
handed metamaterials, and electron nanodevices with se- 
lective transmission and/or reflection. 

F.M.I, acknowledges partial support by the VIEP 
BUAP-2008 grant. 
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